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Hydrodynamics of Cathode Films 

Garbis H. Keulegan 

A hydrodynamical theory of the cathode diffusion layer based on the assumption of 
a uniform density reduction at the cathode surface is presented. By means of an approxi- 
mate solution for the flow and the diffusion equations, the layer dimensions, the velocity 
distributions, and the densities are determined. Further derivations from the solutions 
directly applicable to results obtained by the freezing method are effected. The data of 
Brenner on copper and nickel baths are considered. Using the observed values of the 
density reductions at the cathode surface, the maximum width of the layer, the average 
densities for a given distance from the cathode, and the mobilities of the cations are deter- 
mined and compared with the observations. The agreement between observed and com- 
puted values is regarded as fairly satisfactory. 



1. Introduction 

Neither the cognizance of the phenomenon repre- 
senting the subject matter of the present paper, nor 
the statement of the basic causes producing it is 
original with us. Some 12 years ago, W. Blum, 
Chief of the Electrodeposition Section, National 
Bureau of Standards, described to us the cathode 
film phenomenon as a thin layer of moving material 
very close to the cathode surface and as being 
caused by the continued reduction of the densities 
due to electrodeposition. Since the phenomenon in 
question is also partly hydrodynamical in nature, 
he suggested that we should prepare a theoretical 
study with a view to determining the density dis- 
tributions in cathode films. 

Blum's first thought was that the film might con- 
stitute a seat of discontinuity in densities, since to 
the naked eye a line of demarcation between the 
contiguous layers of liquid near the cathode surface 
is distinctly visible. The line of demarcation 
should be due to the different indices of refractions 
of the liquid layers on the two sides of the line. Thus, 
according to these ideas the flow near the cathode 
surface is in one or two layers having dissimilar 
densities with sharp discontinuities. But, despite 
the appearance of the region close to the cathode 
surface, the opinion was ventured on our part that 
the exsitence of the sharp discontinuities could not 
be shown on the basis of purely hydrodynamical 
reasoning. The basis for this opinion was the belief 
that the phenomenon in question was essentially 
similar to the type of flow encountered when a 
heated vertical plate in the presence of cold air 
sets up convection currents. Since experience in 
this particular case does not show any discontinuities 
in temperatures, it might be expected that density 
discontinuities would not be present in the cathode 
plate problem. Unfortunately, due to other pressing 
work at that time, it was impossible to proceed 
with a mathematical solution of this interesting 
problem. 

About 3 years ago, Blum repeated his request for 
an analytical study of the problem, since he was 
planning to resume the experimental investigation of 



cathode films. The computations that appear in 
this paper were completed shortly after the receipt 
of this request. Meanwhile two other papers on 
the subject, one by Agar fl] ^ and the other by 
Wagner [2], were published. Since the present 
work proceeds on a mathematical basis different 
from these two papers and since the work was 
conceived with a view of interpreting the results 
of A. Brenner, a member of the staff of the Bureau, 
on cathode film densities, it is hoped that some interest 
in the results achieved will develop. 

In the present treatment it is assumed that the 
reduction of densities over the vertical cathode 
surface is constant. On the whole the treatment is 
self-contained and is given in full detail. In general, 
the results shown are approximate and sufficient, 
it is believed, to evaluate average densities in a 
vertical plane parallel to the cathode surface. In a 
future publication it is proposed to investigate the 
problem of cathode films with the basic assumption 
that the density gradient at the cathode surface is 
constant and uniform, rather than the density. 

2. Molecular Character of Electric 
Conduction in Solutions 

Certain fundamental relations of electrochemical 
import that will be considered during the course of 
the analysis are based on the ideas of Larmor [3] in 
the matter of electrolytic conduction. 

Taking an idealized situation it will be supposed 
that the electrolytic solution is a binary system and 
that ionization is complete. In each portion of the 
solution the number of the positive ions and the 
negative ions is the same and the conditions are 
steady, that is, independent of time. At a point of 
the solution let the number of ions of each kind 
per unit volume be n. Measuring positive y horizon- 
tally away from the plane cathode surface toward 
the anode, let a section normal to y be selected. 
Let dNi be the number of the cations and dN2 the 
number of anions crossing the section per unit area 
during time dt. The motions of the ions are due in 

1 Figures in brackets indicate the literature references at the end of this paper. 
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part to the electrical force dP/dy, P being the elec- 
trical potential, and in part to tlie forces producing 
diffusion. If Ui be the specific mobility of the cations 
and ki the coefficient of diffusion of these ions, 



^'K'^^'f+^'S)'^^' 



and for the anions, in similar manner, 



(1) 



(2) 



Since there is no accumulation of electricity at a 
point, the relation between the electrical current 
intensity, I, and the motion of the ions will be 



se{dNi + dN2) = Idt, 



(3) 



where se is the charge on each ion. Applying the 
principle of continuity of the ionic flows, 



1 d fdN, 



dN, 



2 d7j V dt dt 



y 



dn 
~di' 



(4) 



On the other hand, the current intensity remaining 
constant, from eq 3 



d'Ni , (iWi 



0, 



dydt dydt 
and this together with eq 4 gives 



and 



d^Ni 


dn 
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(6) 
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Differeiitiatino; (H) 1 and 2 with respect to y and 
using the relations just, given, 



and 



dn 
Hi 



jj d { dP\,, dh 



d^n 

V 



dn 
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i=-^4("f)+^w^- 



(8) 



(9) 



Multiply the first of these by 1j2, the second by Ui, 
and adding the resulting equations, one obtains 



dn_k2Ui + kiU2 d^n 
~di~ U,+ U2 df' 



(10) 



Now the form of this equation is exactly the form of 
the diffusion equation of an ordinary nonelectrolytic 
solution. Hence, the changes in the concentrations 
of an electrolyte proceed with the diffusitivity con- 
^ stant 

Ul+U2 ' ^ ^ 




Figure 1. Geometrical and physical quantities relating to the 
diffusing layer. 

As this relation holds good however slight the electric 
current may be, it may be supposed that the relation 
holds in the limit when the current is nil. Therefore, 
D is the ordinary coc^fficient of diffusion of the solu- 
tion, and thus a relation exists coimecting this coeffi- 
cient with tlie ionic quantities, that is, with the spe- 
cific mobilities and the coefficients of diffusion of the 
two ions. 

If the relations 



k, = CU, and k2^CU2 



(12) 



exist, where C is a constant independent of the con- 
centration of the ions, then 



D= 



2U,U2 fj 



Ui+U, 



(13) 



Now, from the theory of electi'olytic diffusion by 
Nernst [4] we have 



^~Ur+U2 ' 



(14) 



where T is the absolute temperature and 7? the gas 
constant. Accordingly, the constant C in eq 13 
may be identified with the product RT in eq 14. 

3. The Equations of Motion and of Diffusion 

Taking the case of a plane cathode surface, held 
vertically, a point at the lower edge may be chosen 
as the origin of the rectangular coordinate system 
{x,y) (see fig. 1). The positive x axis is drawn 
vertically upward and the positive y axis normal to 
the surface and into the electrolyte body. Denoting 
the velocity components along the x and the y axes 
by u and v, respectively, the equations of motion are 



c)u , du 

u^ — h^^-= 
ox oy 



1 bp 
p dx 






(15) 



and 
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0--- 



1 d^ 



(16) 



where p is the pressure, p the kinematic viscosity 
of the sohition, p the density, and g the constant of 
gravitational acceleration. These are the equations 
of motion used in viscous boundary layer theory. 

The pressure must be expressed in terms of the 
densities. Denoting the density of the electrolytic 
bath at points far removed from the cathode surface 
by pi, the density near the surface may be expressed 
by the relation 

p-Pi-Ap, (17) 

where Ap is the deficiency of the density in the layer. 
The greatest reduction of density is at the cathode 
surface, that is at the plane ^ = 0. Let the average 
value of Ap over the entire length of the cathode sur- 
face be denoted by Ap^. The study of the densities 
in the la^^er is done conveniently, as will be seen 
later, by considering the ratio 



^ Ap 
Apm 



(18) 



Another quantity important for the analysis is the 
effective gravity 

g^^=G, (19) 



pi 



Now, the second boundary layer equation, eq 16, 
states that in the layer next to the cathode surface 
the pressure along a horizontal plane remains con- 
stant. Moreover, all motions being nil in the main 
body of the bath solution, the pressures are hydro- 
static and hence 

p=Po-pigx, (20) 

where po is the pressure at points in the horizontal 
plane passing through the origin. Thus, 



1 ^__ Pi. 
p dx p 



(21) 



Remembering that Ap is a small quantity and mak- 
ing use of the relations in eq 18 and 19, we have 



p dx 
and thus eq 15 reduces to 



u^ — [-V -^—= 
ox oy 



■g = Ga, 



dy 



, + Ga. 



(22) 



(23) 



This is the equation of motion that is appropriate to 
the problem at hand. To this must be added the 
conditions of continuity, 



dx dy 



(24) 



Next will be considered the equation of diffusion. 
In the presence of convective currents, the diffusion 
equation takes the form 



bn , bn , dn -^/ 






(25) 



where n is the number of particles of the diffusing 
material per unit volume. This relation reduces to 
Fick's equation when convective currents are absent. 
Although the equation is well known, it may be help- 
ful to give a derivation in order to clarify the ideas 
involved. Neither the name of the investigator 
giving the first derivation nor the manner of deriva- 
tion is known to us. 

Consider a parallelepiped of sides dx, dy, and 1. 
Let the number of particles crossing the face at x per- 
pendicular to the X axis during time dt be 

Nudydt, Nu = -D^+nu, 

and the face B^t x-{-dx be 

Nudydt=fNn+^^dx^dydt. 

Similarly, let the number of particles crossing the 
face at y perpendicular to the y axis be 



N21 dx dt, 
and the face 2iiy-\-dy be 
N22dydt=(N2i 



■AT 7->^^ . 



biV21 



by 



dy jdxdt, 



so the accumulation of the particles in the parallele- 
piped is 



dx dy 
But the accumulation is also 
dn 



dt 



dx dy dt. 



Hence, 



dn_ dNu dN2i 



dt 



dx dy 



or 



dn 
~dt 



-( 









dn d7i 

-u^ z;— - 

ox dy 



-'< 



du dv\ 
d^^d^J' 



The last term in the right-hand number vanishes 
because of the incompressibility of the solution. 
This now proves the general equation of diffusion 
given b}^ eq 25. A tacit assumption of the derivation 
is that the ratio of the volume of the diffusive 
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particles to the volume of the solution is a small 
fraction. Stated otherwise, the result is valid for 
small concentrations. 

The diffusion equation is m,ucli simplified when 
applied to the cathode film, boundary layer. First, 
the concentrations do not vary appreciably with 
time. Second, the layer being very thin and com- 
paratively very long, the quantity d^n/dx^ is sm,all in 
com.parison with b^n/dy^. Accordingly, eq 25 becomes 

bn c>n_ ^ b'^n 
dx by dy^^ 

and since n is proportional to a, also 



da , i>a 

u-:^ — h^ ^-= 

ox oy 



D 



by' 



(26) 



This is the diffusion equation appropriate to the 
problem at hand. Equations 2:^, 24, and 26 govern 
the m.ovem.ent of liquid and i\w distribution of the 
densities in the cathode lihu subject to the boundary 
conditions. 

4. Boundary Conditions 

The boundary conditions to be considered are of 
two kinds. One kind relates to the velocities and the 
other to the densities, the latter having electro- 
chemical implications. 

The boundary conditions to be assigned to vc^locity 
components are the most obvious ones. The solution 
or the liquid of the film being viscous, u and v vanish 
at the cathode surface ; t hat is 



u=v=0, y = 0. 



(27) 



To state the condition for the velocities at large dis- 
tances away from the cathode surface it is necessary 
to introduce first the concept of the boundary layer 
thickness (see fig. 1). The boundary \n,yeT thickness, 
ysj is such that the density of the solution at this 
distance from, the cathode differs only im.perceptibly 
from, the density of the bath. Since the m.ovem.ent 
of liquid in general is nil in the m.ain portion of the 
bath, a second condition for the velocities is 



u=o. 



y=ys' 



(28) 



It is obvious that the thickness t/s of the boundary 
layer at the lower edge of the cathode surface, that is 
at .1=^0, vanishes. This m.eans that the liquid m,oving 
upward in the film, m.ust be replenished by a later 
and horizontal movement of the bath solution into 
the film, region. Hence 



v>0, 



y=ys' 



(29) 



The assum,ption of vanishing velocities in the bath 
in the vertical direction is an idealization. Under the 
actual conditions of experimentation, arbitrary m,ove- 
ments in the bath may be present even when forced 



agitation is not resorted to. Since the m.ovem.ent of 
liquid in the film next to the cathode surface is up- 
waid and in the film, next to the anode downward, 
it is obvious that a circulation is induced in the bath 
liquid. Due to accidental variations of temperature 
at places the regularity of the circulatory motion will 
be disturbed and will give rise to weak eddying mo- 
tions. These will be most noticeable^ near the upper 



surface of the bath and close to the cathode 



5uch 



accidental disturbances can be imagined to have only 
m.inor effects on the main motion in the film, and thus 
will be ignored. 

The boundary conditions that are to be assigned 
to the densities are not obvious, especially the one 
relating to the value of the densities at the cathode 
surface. Regarding this, it seems that two different 
statements may be made depending on whether the 
transmitted electric current is of the limiting value 
or is inferior to this limiting value. Taking the case 
of weak electric currents, we may considcu* the ques- 
tion of the densities at the cathode surface as follows: 

Assume that the ionic system is binary and 
ionization is complete. Since at the cathode surface, 
that is at ^==0, dN2/dt=0, eq 3 gives 

dN, I 



dt se 
and thus eq 1 and 2 now become 



and 



I ^^ dP , J dn 

se dy ay 



^ ^j dP . dn 
^^''^^dy-^^dy' 



where n is the number of ions of each kind per unit 
volume. Eliminating the potential gradient between 
these two equations and making use of eq 11, 



n ^^_ 



C/2 



dy Ui+U2se 



(30) 



This is the relation between the electric current 
intensity and the gradient of n at the cathode surface. 
It is desirable that the quantity n in eq 30 be 
replaced by c, the latter being defined as the gram 
equivalent of the salt per cubic centimeter of the 
solution. In a unit cube there are n cations and n 
anions each of ionic valence s, and the system is 
completely ionized. Introducing N, the Avogadro 
number, the molar concentration is n/N and hence 
the concentration in gram equivalent of the salt per 
cubic centimeter of the solution is 



ns 



and eq 30 now becomes 



D 



dc 



U, 



dy Ui+U, eN 



(31) 



957570—51- 
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Since eN is the universal constant F, with the 
numerical value F= 96450, and with the dimensions 
coulomb per gram equivalent, 



V 



dc 



U, 



dy U. + U^F 



(32) 



The same condition may be expressed in terms of 
densities instead of concentrations. For the major- 
ity of electrolyte solutions the variation of densities 
with the concentration is linear; that is 



dc 



-=m. 



(33) 



For copper sulphate m= 1.32X10"^ gram equivalent 
per gram; and for nickel chloride m= 1.71X10"^ 
gram equivalent per gram. Equation 32 may now 
be written as 



mD 



U, 



dy-U, + U,F' 



y = 0. 



(34) 



This is the relation between the intensity of the 
electric current and the gradient of the density at a 
point of the cathode surface. It is the only relation 
that may be derived on the basis of simple concepts 
of an electrochemical nature. Accordingly, we are 
not certain of the state of the densities at the cathode 
surface for the present. 

Two simplifying assumptions may be made. First, 
it may be supposed that the current intensity is 
uniform at the cathode surface. The resulting 
boundary condition would be that the density gra- 
dient is uniform at the cathode surface. Second, it 
may be supposed that the density is constant at the 
cathode surface. For the present analysis we shall 
adopt the second alternative simply for the reason 
that the corresponding analysis may be completed 
with ease. Thus it will be supposed that, 



and 

In view of eq 18 

and 



Ap=0, y=^. 
a = 0, y= oo . 



(35) 
(36) 



The set of equations 27, 28, 35, and 36 constitute the 
boundary conditions to be considered in solving the 
set of equations 23, 24, and 26. 

Since the relation implied in eq 34 is valid for all 
conditions, it may be used to consider average effects. 
Expressing p in terms of a through the use of eq 18, 
multiplying the two sides of the resulting equation 
by dx, and integrating between x^O and x = /i, where 
Zi is the length of immersion of the cathode surface, 



-mPAp, 



ni /da\ 
Jo \dyj 



dx = 



U, 



U,+ U, F' 



(37) 



where now / is the average intensity of the electric 
current measured in amperes per square centimeter. 
We shall use this relation to compute the transport 
number of copper ions taking our data from the results 
of the Brenner experiments. 

5. The Dimensional Form of the 
Basic Equations 

The solution of the basic equations subject to the 
boundary conditions indicatecl above will be greatly 
simplified if the equations are changed into forms in- 
volving dimensionless quantities only. The clue to 
the proper course of analysis may be found in the 
method of Blasius for treating the viscous boundary 
layer of plates in a current of air. 

The limiting values of the densities and the veloc- 
ities being uniform at the boundaries of the cathode 
film or of the diffusion layer, it is to be expected that 
the distribution of the velocities and of the density 
reductions in different normal sections are affine to 
each other. The affinities will be established by the 
introduction of the characteristic length 5 and the 
characteristic velocity [/, both of which are functions 
of X only. 

Introducing the stream function \l/, such that 



u= 



dy 






(38) 



the condition of incompressibility, eq 24, is automa- 
tically satisfied. In terms of the stream function the 
equation of motion, eq 23, and the equation of diffu- 
sion, eq 26, transform to 



and 



d^ly dV ^^ ^V c)V , ^ 



dy dydx dx dy^ 
d\l/ da dxp da 



D 



d^c. 



dy dx dx dy dy^ 

Introducing the new variables 



and 



H- 



m' 



(39) 
(40) 

(41) 
(42) 



and remembering that U and 8 are functions of x 
only, we have, first 



and 



and hence 



d__dri d_\d_ 
dy dy drj 6 dr]' 

d _dT] d _ 7) d8 d 
dx dx dr] 5 dx d-q' 



dx dx dx dr)' 



(43) 
(44) 

(45) 
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M_^ dH 
by dr,' 


(46) 


dV dUdH UvdBdm 
)a; by dx dr, 5 dx (iij^ 


(47) 


dV U d'H 

by' 8 dr,'' 


(48) 


dV udm 

by'' 5' dv' 


(49) 



and 



Substituting these expressions in eq 39 and 40, 

dx\dr, / 5 dx ' dr,^ 5' dr,'' ' 

(50) 

(51) 
(52) 



and 
where 



d dx^ ^ dr] ' b^ dr}^ 



V 

D' 



Now eq 50 and 51 will be freed of the differential 
operators involving x if one can write 



and 



dx 8 ax yb^ 



d dx^ ' 7^ 



-G 



(53) 



(54) 



where Mi and M-i are numerical constants. In view 
of tliese provisional lolations, the two equations 
simplify to 



dm 



and 



dn-" 



.Ijrd'H 1 (dHV, . ,.,, 



Mi dr, 
H da 



d^l 

M2dr,'^dr,'~~ 



-0. 



(56) 



In fact, the provisional equations 53 and 54, 
determine the dependence of the characteristic 
velocity U and the characteristic length 5 on the 
variable x. The provisional equations 53 and 54 
imply 



^^^S=^^^(^^)' 



and 



d 

{Ud) = 
7 

vU=yb'G. 



M^bf^m 



(57) 



(58) 



(59) 



It is now obvious that U and b each are single 
power terms of x. Assuming 



U^Ax\ 



(60) 
(61) 



where A and B are constants involving tlu^ physical 
characteristics of the electrolytic solution. Substi- 
tuting these expressions in eq 57, 58, and 59, com- 
paring the powers of x and the coefficients for the 
individual equations, the com.parisons yield 



and 



M,r=M2(s-\-r), 

M2(s + r)AB'=v/y 

2s+r—l = 0, 

r=2s, 

vA=yB'G. 



Equations 64 and 65, yield r=l/2 aiul 
Equation 62 now yields 

2ilfi— 3M2=0. 



(62) 
(63) 
(64) 
^65) 
(66) 

-1/4. 



For sim,plicity subject to the last derived o(i nations, 

we may put Mi =1/2 and il/2=l/3. 

Substituting the values of the constants A and B 
that are o})taine(l from iH\ 63 and ()6 into eq 60 and 
61, there results 

U = 2G'/'x'^', (67) 

and 

/4j,2\1/4 

which are the functions showing the dependence of 
the characteristic velocity U and the characteristic 
length b on the independent variable x. 

Substituting the values of Mi and M2 as found 
above into eq 55 and 56, these simplify into 



^^.?H-f^- 



and 



■Kf)"+«-. 



5+3H^-0, 



dr 



(69) 



(70) 



which are the equation of m.otion and the ecpiation 
of diffusion expressed in dim.ensionless forms. 

In term.s of the new variables, the velocity com- 
ponents u and V are, using eq 38, 45, and 46, 



u= 



and 



drj 



■'-Tx^'''^''-^d-x- 



dS dH 



dr, ' 



(71) 
(72) 
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and the boundary conditions are 



dH 

dr] 



= 0, V-- 



(73) 






The two equations 69 and 70 becon),e the basic 
equations for the problem, of tem.perature variation 
of a heated plate placed vertically in a cool air bath, 
when 



_T-T\ 



(74) 



where T is the temperature near the plate, To the 
tem.perature of the plate, and Ti that of the cool air, 
and the parameter 7 is replaced by 



r=vlk, 



(75) 



where k is the thermometric conductivity. 

Goldstein [5] has discussed the Pohlhausen solution 
of the heated plate problem, and the equations con- 
sidered are similar to eqs 69 and 70, when allowances 
are made for the variations in some of the coeffi- 
cients. Thus, in the two basically different phe- 
nom^ena of the diffusion layer near a cathode plate 
and of the conduction layer near a heated plate, 
the corresponding parameters 7 and <t play identical 
roles. The latter is well known as the Prandtl 
number. The form.er will be called the diffusion 
param.eter for the present, since it is a dim.ensionless 
number expressing the ratio of the kinem.atic vis- 
cosity to the coefficient of diffusion. In forming 
the numerical value of 7, both D and v m.ust be 
expressed in square centim.eters per second. 



6. An Approximate Solution of the Basic 
Equations 

Since an approximate solution is adequate for 
the present objectives, it is hardly necessary to 
effect a rigorous solution of the two equations, 69 
and 70. To prepare the way for the approximate 
solution, it is necessary to introduce the boundary 
layer thickness ?/s, which is a function of x. To 
the distances y and ys correspond the parametric 
distances 7; and tj^, such that 



T]=ylb and Vs^Vs/d, 
We introduce the ratio 

o=vlvs=ylys' 



(76) 



(77) 



It is obvious that the quantity d will serve as the in- 
dependent variable in the place of 7]. Indeed, the 
basic equations 69 and 70 transform to 



7^+3,.H^^-2,J^j+,3,^0, (78) 



and 



-2+377, /i -J-^^O, 



de' ' "^'"^^ dd 

with the boundary conditions 

H=0, 9=0, 

6=0, 



(79) 



dH 



de ^' 



1, 



(80) 



a=0, 6=1, 

in place of the set shown by eq 73. 



It is obvious, intuitively, that the space rate of 
change of the velocity components, bujby a.nd dv/dyy 
each vanish at large distances from the cathode sur- 
face. These conditions are fulfilled likewise quite 
approximately at distances equaling or exceeding yg. 
This observation suggests the additional boundary 
conditions 



d'H 
dd^ 

da 
dd 



= 0, ^=1, 



= 0, ^=1, 



(81) 



the first of these resulting from eq 48. All the 
boundary conditions of eqs 80 and 81 are satisfied 
if it is assumed that 



a= 1 —sin 



2 ' 



and 



H- 



A, 



4.1 r 



-cos -kB—— cos 27n 
4 



a 



(82) 



(83) 



The above analysis involves the quantities Ai 
and ris, which at present are not known. For the 
determination of their values the differential equa- 
tions 78 and 79 are at our disposal. Multiplying 
these equations by dd, integrating between the 
limits ^=0 and ^=1, and making use of the boun- 
dary conditions, the result is 
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and 



(a+-/: 






(85) 



after resorting to the rules of integration by parts. 
Introducing the values of a and H from eq 82 and 
83 and effecting the integrations, we find 



-2x7A- 



25 



i^H?(^)^ 



= 0, 



T, 14 Ai 



and hence 
and 



4 5^* / ,125 \ 

'''-14(x-2)V^ + 448 V 



A- 



5x2 



287,, 



(86) 
(87) 



The expressions for r^^ and Ai as givc^n by eq 86 
and 87, respectively, and the expressions for a 
and H as giv(^n l)y eq (S2 and 83, respect iv(^ly, con- 
stitute the first approximate^ solution of the (liffusioii 
layer problem based on the assumption that the 
reduction of densitic^s, or concentrations, at the 
cathode surface is imifoim. 



7. Thickness of the Cathode Diffusion Layer 

Inasmuch as the difi'usion parameter y=v/D is or- 
dinarily a ver}^ large quantity, it is sufficient to write 
in the place of eq 86 



Stt^T 



14(x-2) 



(88) 



Introducing the value of ys from eq 76 and the value 
of d from eq 68, wc find first 



Vs 



r 5x^7 T^T4z^n^/^ 

Il4(7r-2)J [ ~ 



LGyU 



and next introducing the value of G from eq 19 



:3.33 



\gApmyJ 



^1/4 



(89) 



This is the expression that gives the thickness 
of the diffusion layer along the cathode surface 
(See fig. 2). The thickness is not uniform., but 
increases from, the zero value at the lower edge of 
the cathode surface to a maximum, at the upper end 
of the cathode surface still in contact with the 
electrolyte. The rate of increase of the thickness, 
however, is very gradual due to the fact that the 
thickness if proportional to the one fourth power 
of X, the distance from the lower end of the cathode. 
Again at a chosen point x the thickness of the layer 
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Figure 2. Shape of the boundary of the diffusion layer. 

decreases when Ap.^^, the m.axim.um. reduction of the 
density at the cathode surface, is increased. Since 
the decrease of density at the cathode surface is 
proportional to the electrical current (hmsity, when 
the latter is inferior to its limiting value, it is to 
be expected that a decrease in the intensity of 
current will cause the thickiu^ss of the difFusion layer 
to be augment (m1. 

Having established the* magnitude of the (hfl'iision 
layer thickness, it is of interest to inquire next into 
the nornuil density gradient at the cathode surface. 
The variation of the density gradient will be best 
understood if it is expressed in terms of the average 
density gradient, the averaging being made over the 
entire cathode surface. In view of eq 18, the quan- 
tity 



KbyjJKby), 



11 



(90) 



where the luuneiator on the left-hand side refers to 
local values at the point x and the denominator to 
the average values, is the ratio to be considered. 
Since 

da_c>a d£ 

c>y be dy^ 



we have from eq 77 and 82 



TT 



or 



where 



Ci- 



K^yJo x" 

6.66\^Ap^ 



-1/4 



(91) 



(92) 



Accordingly, the average value over the length li 
of the cathode surface is 



( 






(93) 
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and the desired ratio is 



B- 



©"' 



(94) 



which gives the proportionate variation of the normal 
density gradient over tlio cathode surface The 
variation is shown graphically in figure 3. In the 
section of the cathode from x/li = 0.2 to x/li = l.O the 
ratio is nearly constant. This circumstance assures 
the utility of the present solutions even for the case 
in which the normal concentration gradient is strictly 
constant over the entire surface of the cathode, pro- 
vided that one deals with the average values of the 
densities in the planes close and parallel to the 
cathode surface. 

If the diffusion layer thickness at x=li, that is at 
the upper end of the cathode surface just in contact 
with the hath, be denoted by ysi, then the average 
value of the normal density gradient, expressed in 
terms of the maximum density reduction, is 



(95) 



/da\ _ 2.1 

\dyJo~ Vsi ' 

Use will be made, subsequently, of this relation. 
8. Velocitus in the Diffusion Layer 



Consider first the velocity component u, that is 
the component parallel to the surface of the cathode. 
Now, 



where 
Since 



U=LJ -j-y 

a 7] 

dH^dH dd 

dr] do df] 



(71) 
(67) 
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Variation of the density gradient along the cathode 
surface. 



eq 71 may serve, in view of eq 77, 83, and 87, to 
give the distribution of the velocities in a normal 
section of the diffusion layer. There is, however, 
another method for representing the velocity dis- 
tribution. Moving away from the cathode surface 
in a normal section, u is zero at the cathode surface, 
increases to a maximum value Um, and then de- 
creases to zero at the limit of the boundary layer, 
that is at y=ys- Thus, the distribution of the 
velocity may be obtained also by considering the 
ratio u/Um- Now 



"-=^(^1- 



In terms of the variable 6, 



/dH\ ^}_fdH\ 

\dri /„ Vs\dd/m 



(96) 



(97) 



If ^1 be the value of 6 that makes u a maximum, we 
have from eq 83, 



/dH\ 



-Ai sin 7r^i+-sin 27r^i |- 



The value of 6i = % satifies the relation. 



cos TT^i+COS 27r^i==0. 



This makes 



sin 7r^i+- sin 27r^i = 



a. 299, 



and therefore 



Kdd/rr 



a. 299 A 



(98) 



Substituting into eq 96, from, eq 98, 97, 88, and 67 
gives 



'^.. = 0.821 



or using eq 19, 



<f) '■"■ 

^gx Vpi7/ 



(99) 



Accordingly, the m.axim,um value of u in a normal 
section changes from, one section to another and 
increases with the one-half power of the distance 
from, the lower end of the cathode. 
Since, 

u/Um. = {dHld7])l(dHldrj)m . 

and introducing the relations in the equations 97 and 
98, the distribution of the velocity component u in sl 
normal section in the diffusion layer becomes 

u/Um = 0.77 TsinTT^+^sin 2t^\' (100) 
L Vs 2 ysj 
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Figure 4. Distribution of the vertical velocity component. 

The graphical representation of the distribution is 
given in figure 4. 

Next consider v, the vekx'ity component noi-iual lo 
the cathode surface. The velocit}^ can be written 
from eq 72 as 

The m,agnitude of v is hirgest at the liquid boundary 
of the diffusion layer. Let Vm be this largest value. 
This occurs at tj^Vs. Let Hg be the value of // at 
^ = 18 or at 6=1. At t] = 't)s^ dll/(lr]=0, and thus 



where 



dx ' 



(102) 



Hs = 2A,iTr. (103) 

In view of the values of U and 5 from, eq 67 and 68, 



f7g=l(4.^^)-. 



-1/4^-1/2^ 



Substituting the above result in eq 102 and also 
introducing the value of Ai from, eq 87, 



,= -1.013(^2^)^/^7"'^'^"'^'. 



(104) 



Accordingly, the liquid of the bath moves hori- 
zontally into the cathode diffusion layer. The veloc- 
ity with which the liquid approaches the cathode 
surface is variable along the length of the cathode 
surface, the m,axim,um. value occurring at the lower 
end of the cathode surface. After a certain elevation 
is reached the variation in the value of Vm ceases 
to be pronounced. 
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FicaRio 5. Distribution, of the horizontal velocity component. 

The representation of the distribution of v in a 
normal section of the layer is most conveniently 
doiu^ by considering the ratio i^lvm. Now, 

or in tc^rnis of ^, 

Vm H.L 3 de_\ 

This gives, in view of eq 83 and 103, 



-- cos ird—— cos 27r6 
2 o 



■ - sm TT^ — y- sm 27r^. 
(105) 



The graphical representation of the distribution is 
shown in figure 5. 

9. The Brenner Freezing Method and the 
Laminar Mean Densities 

In the Brenner method [6] the electrolytic solution 
near a cyhndrical cathode is suddenly frozen, succes- 
sive layers of the solid of 0.003- to 0. 004-inch thick- 
ness are carved off on a lathe, and the compositions 
of the layers are determined. Thus, in this method 
the average concentration or the average density of 
a lamina of thickness diji, of distance i/i from the 
cathode and of length /i are determined. Accord- 
ingly, if theoretical results are to be compared with 
the observed results of the Brenner experiments, it 
becomes necessary to derive a formula giving the 
average value of densities of a lamina of distance 
yi from the cathode surface, 
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Figure 6. Geometrical and physical quantities relating to the 
diffusion layer in Brenner^ s tests. 

Let the average density reduction, in the lamina of 
length Zi and of distance i/i from the cathode surface 
be denoted b\^ Ap\ Mathematically, 



1 r^i 

Ap^=Y Ap dx. 



(106) 



Let Iq be the value of x for which the boundary layer 
thickness equals t/i (see fig. 6) . Obviously, for values 
of x<Clo, Ap = Q, and thus, 



^'Hi'!^"^''- 



(107) 



Since, according to the present solutions, the density 
reduction at y=yi is given by the relation 

Ap = Ap„(l-sin|^) (108) 

(from eq 82), then 

Ap'=Ap.[;-J^.-y\inf|;..]. (109) 

It must be remembered in effecting the integration 
of the second integral of the right-hand member, that 
2/1 is held constant. Putting 



and 



yi/ysi=Pi, 



(110) 



It is seen from the law of the layer thickness, eq 89, 
that 



and thus. 



x = lS\ 
dx = —4:lodp/p^. 



(Ill) 
(112) 



(113) 



Accordingly, in terms of the new variable fi, 

Ap'=Ap^[l-^t + 4/3tJ^''sin|/3^. (114) 



For convenience, we may write 



F{p,) = P\-4p',j\m^^ 



dp 



and hence. 



Ap'/Ap^=^-Fip,). 



Since Pi = yilysi, we have also 



Ap7Ap^-l-i<^(7/i/2/.i). 



(115) 
(116) 

(117) 



This is the desired formula that will be used in 
analyzing the experimental data of Brenner. It 
involves the average reduction of density over a 
lamina, the average maximum reduction of density 
at the cathode surface, the distance of the lamina 
from the cathode, and the maximum thickness of 
the diffusion layer. The determination of the 
function F(Pi) was carried out by replacing the term 
sin (7r/2)/3 appearing in the integrand by a power 
series in 13 and then effecting the integration. The 
results of this determination are given in table 1. 

10. An Analysis of the Results Obtained by 
Brenner with the Freezing Method 

The comparison between the results of the present 
analytical theory of the cathode diffusion layer and 
the results of the Brenner experiments on the 
reduction of the densities in the diffusion layer will 
be carried out along the following lines. First, from 
the observed values of the reduction of density at 
the cathode surface, Ap^, the maximum thickness 
of the diffusion layer, ysi, will be computed and 
compared with the observations. Second, using the 
same observed values of Apm, the ration ApjApm 
will be computed for the distance yi, and these will be 
compared with the observations. Third, Ap^ will be 
computed from a consideration of the electric 
current intensities used in the expeiments and will 
be compared with the observations. 

Table 1. Numerical values of the function F (yilysi) 



yilvsi 


Fiyr/Vsi) 


y\lvx\ 


Fiyifysi) 





0. 0000 


0.6 


0. 8926 


0.1 


.2070 


.7 


.9388 


.2 


.4024 


.8 


.9640 


.3 


.5708 


.9 


. 9796 


.4 


.7099 


1.0 


1. 0000 


.5 


.8180 







The Brenner [6] results relate to electrolytic baths 
of copper sulfate and nickel sulfate. The baths 
actually used were not pure binary systems, but 
this fact will be ignored. The complete data from 
the Brenner tests are collected in table 2. The 
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Table 2. Diffusion layer densities of the Brenner test 



Distance 
from 


Normal- 
ity, 


Density, 


Current 
intensity, 


Metallic 


cathode, yi 


N 


P 


/ 






Gram 








cm 


equivlliter 


«7/cw3 


nmp/cm^ 




3.69X10-2 


1.274 


1. 1009 






11.16 


1. 522 


1.1174 






18.79 


1.800 


1. 1381 


[ 0.04 


Copper 


27.63 


1.920 


1. 1465 






37.63 


1.931 


1. 1473 






3.83 


1.741 


1.132 


1 




11.18 


1.885 


1.142 






21.05 


1.971 


1.149 


[ O.Ol 


Nickel 


32.46 


1.985 


1.150 






44.12 


1.985 


1.150 






3.83 


1.561 


1.118 






12.36 


1.724 


1.130 






21.05 


1.872 


1.139 


[ 0.02 


Nickel 


30.00 


1. 935 


1.146 






40.00 


1. 959 


1.148 






3.83 


1.382 


1. 105 






12.36 


1.694 


1.128 






21.05 


1.900 


1. 143 


I 0.04 


Nickel 


30.00 


1.947 


1.147 






40.00 


1.981 


1. 149 







concentrations in the original data are in terms of 
the normalities. Tlie dc^nsities as (Ieriv(Hl from llie 
normality values are added. The avei'age densities 
p' taken from the table are plotted against iji in 
figure 7, and smooth curves are drawn through the 
points of observation. The extreme values poand pi 
of the densities are read from the curves. The 
differences (pi — Po) give Ap,n, the average* maximum 
reduction of densities at the cathode surfaces. The 
difference (pi — p'), where p' is the observed value 
of the densiti(?s corrt^sponding to the distance ?/i 
gives Ap\ the average reduction of the (k^nsities 
corresponding to the distance ?/i All these quan- 
tities, together w^itli the ratio Ap'/Ap,^, are entered 
in table 3. 

From the observed values of Ap^ and pi the nnixi- 
m.um. layer thickness, y^i, m.ay be computed from 
the theoretical formula in eq 89. However, the use 



of eq 89 for this purpose requires that the parametric 
diffusion number, y, be known. This number in- 
volves the ratio of the kinematic viscosity p, to the 
diffusion constant D, for each of the two electrolytic 
solutions. The viscosities of the electrolytes may 
be determined from the normalities using the 
Arrhenius [7] formula, 

m=pA 

where pi is the specific viscosity of a normal solution 
referred to the solvent at the same temperature, 
p is the specific viscosity of the solution having the 
normahty N. Wagner [8] gives pi = 1.358 for copper 
sulfate and pj::^ 1.361 for nickel sulfate at the tem- 
perature 25° C. The normalities of the bath solu- 
tions in the Brenner tests being 2 and the densities 
about p= 1.151, the kinematic viscosities of the two 
solutions are practicahy equal and so p= 0.0 144 
cmVsec. Regarding the (hfl'usion coefFicient D, the 
Thovert [9] values for co])per sulfate may be sup- 
posed to apply as well to nickel sulfate. For copper 
sulfate at a temperature of 17° C, A^==1.95, the dif- 
fusion coefficient is i?=2.66X10"'^ cmVsec. Hence 
the value 

7 = 0.572X10^ 

may be used for the two solutions. As regards the 
cathode immersion length /j, we take it to be 10 cm, 
this being the value communicated to tlie author 
orally by Brenner. 

Using the above-mentioned values of the i)ertiuent 
quantities the n)agnitu(k^s of the maxim uni layer 
thickness were computed from eq 85, and these are 
shown in table 4. Also there are given in the table 
the observed values of the same quantity as read 
from the curves in figure 7. At instances equalling 
the maximann layer thicknesses the average densities 
equal the bath densities. There is som.e (hfficulty 





Table 3. Reductions in 


the diffusion 


layer densities of the Brenner tests 




Distance 


Reduction of 


Reduction of 


Theoretical 


Proportionate 
density re- 
duction, 

l-Ap'/Apm 


Proportionate 


Current 


Metallic 


from cath- 
ode, yi 


density at 
cathode, Apm 


density in 
layer, Ap' 


maximum 

layer thickness, 

2/«i 


distance from 
cathode, yi/ys\ 


intensity, 

I 


ion 


cm 


g/c cm^ 


glc cm^ 


cm 






amp/cm"^ 




0.00 




0. 061 




0.00 


0.00 






3.69X10-2 




.047 




.23 


.12 






11.16 
18.79 


0. 061 


.031 
.010 


0. 0302 


.49 

.84 


.37 

.62 


0.04 


Copi)er 


27.63 




.002 




.97 


.91 






37.63 




I .000 




I 1.00 


1.25 






0.00 




.025 




f 0.00 


0.00 






3.83 




.018 




.28 


.10 






11.18 
21.05 


0.025 


.008 
. 001 


0. 0382 


.68 
.96 


.29 
.55 


0.01 


Nickel 


32.46 




. 000 




1.00 


.85 






44.12 




. 000 




I 1.00 


1.15 






0.00 




f . 040 




f 0.00 


0.00 






3.83 




. 030 




.25 


.11 






12.36 
21.05 


0.040 


.018 
. 009 


0. 0347 


.55 

.77 


.35 
.61 


0.02 


Nickel 


30.00 




.002 




.95 


.86 






40.00 




I .000 




i 1.00 


1.15 






0.00 




.059 




f 0.00 


0.00 






3.83 




.044 




.25 


.12 






12.36 
21.05 


0.059 


.021 
.006 


0. 0315 


.64 
.90 


.39 
.65 


0.04 


Nickel 


30.00 




.002 




.97 


.95 






40.00 




I .000 




I 1.00 


1.26 







Or.7570— 51- 
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Figure 7. Graphical representation of the data from Brenner^ s 
tests. 

in identifying the exact magnitude of the maximum 
layer thickness from the curves of figure 7. If this 
element of uncertainty be ignored, the agreem.ent 
between the theoretical and observed values may be 
judged to be fairly satisfactory. 

The second question to be considered is the dis- 
tribution of the mean densities in the layer itself, 
that is the variation of Ap^/Apm with yjysi- The 
observed proportional density reductions for a given 
relative distance are given in table 3. The method 
of determinations is obvious from the contents of 
the various columns. In forming the ratio yi/ysi, 
the theoretical values of ^.,1 are used. Values from 
the table are plotted in figure 8. The curve drawn 
is the theoretical distribution and is the plot of eq 
117. Again the agreement between observation and 
theory may be judged to be fairly satisfactory. 

The third question to be exam.ined is the depend- 
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Figure 8. Distribution of densities of the diffusion layer 
according to theory and observation. 



ence of the average density gradient at the cathode 
surface upon the intensity of the electrical current. 
This question deals with the electrochemical bound- 
ary condition represented by eq 37. Since the 
average gradient according to eq 95 is 



\^y)o 

eq 37 m.ay be changed to 
2.1m DApm = 



2.1 
ysi 



u. 



2/sl 



U^+U^F 



(118) 



This relation may be examined in two ways. If the 
transport number of the cations (nickel and copper) 
are known, the relation may be used in computing 
Apm, the average density reduction at the cathode. 
Or, using the observed value Apm, the transport 
number of the cations may be determined. The 
second method will be used here. 

From, the data in table 3 the ratios Apmlysi are com- 
puted and these are plotted against current density 
in figure 9. The distribution of the points is linear, 
and the equation of the line drawn is 



ApJiysJ) = 0.5XlO\ 



(119) 



Table 4. — Theoretical and observed maximum diffusion layer thickness in the Brenner experiments 



Bath 
density, 

PI 


Cathode 
density, 

PO 


Reduction of 
density at 
cathode, 

Apm 


Cathode 

immersion 

length, 

h 


Observed 

maximum layer 

thickness, 

y>i 


Theoretical 

maximum layer 

thickness, 

y>>i 


Current 
intensity, 

/ 


Metallic 
ion 


g/cm 3 
L148 
1.150 
1.148 
1.149 


g/cm 3 
1.087 
1.125 
1.108 
1.090 


g/cm 3 

0.061 

.025 

.040 

.059 


cm 
10 
10 
10 
10 


cm 
0.032 
.032 
.037 
.034 


cm 
0. 0302 
.0382 
.0347 
.0315 


amp/cm 2 

0.04 

.01 

.02 

.04 


Copper 
Nickel 

Do. 

Do. 
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We also have 



and 



m=l. 32X10-2. 



Substituting these vahies in eq 118, it is seen that 



U, 



U.+ lh 



= 0.356. 



(120) 



Accordingly, the transport of the cations, nickel, is 
tic =0.644, a value very close to the generally ac- 
cepted value determined by other methods. Thus 
in this respect also there is agreem.ent between theory 
and observation. 

Another element of the theory that could be ex- 
amined with considerable signficance is the question 
of the vertical velocities. Unfortunately in the ex- 
perim,ents of Brenner this matter was not touched. 
To give an idea of the order of magnitude of the 
velocities involved, we may compute the maximum 
upward velocity in a normal section of the difTusion 
layer using eq 99. Consider the case of copper sul- 
fate. Corresponding to the current intensity of 
1 = 4X10-2 am,p/(^m2, the density reduction at the 
cathode surface is Ap,„/p i== 5.31X10-2. Suice 1/7 = 
1.74X10"^, the form.ula gives for m.aximum velocities 



Ur, 



gx 



^2.50X10- 



Thus at a point a^ = 10 cm., Um equals 2.48 mm/sec. 

I'he discussion presented above must lead one to 
the following conclusion. The approximate analyti- 
cal theory of the diffusion layer based on the assump- 
tion of a uniform density reduction along the entire 
cathode surface givers fairly satisfactory results 
regarding the m.aximum. width of the layer, the aver- 
age densities in planes parallel to the cathode surface, 
and the transport number of ions. Brenner is now 
engaged in further work on the variation of densi- 
ties in planes parallel to the cathode surface. In- 
formation of this kind m.ay possibly point to defects 
in the present theory. In that event it might be 
necessary to base the theoretical investigation of 
the problem, on the assum.ption that the gradient 
of the density and not the density is uniform, along 
the cathode surface. 



The author expresses his sincere appreciation to 
W. Blum, for the guidance shown in the beginning of 
the investigation, and to A. Brenner for advice giv 
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ing the paper. 
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